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BLUNT- BODY STAGNATION- REGION FLOW 
WITH NONGRAY RADIATION HEAT TRANSFER - 
A SINGULAR PERTURBATION SOLUTION 

By Walter B. Olstad 
Langley Research Center 

SUMMARY 

A singular perturbation solution to the blunt-body stagnation-region flow of an 
inviscid, radiating gas has been obtained by means of the Poincare-Lighthill-Kuo, or 
perturbation-of- coordinates, method. A number of results for a gray gas have been pre- 
sented in order to provide some physical insight into the effects of various parameters on 
the shock- layer enthalpy profiles and the radiant heat-transfer rates. 

A nongray absorption-coefficient model was developed which includes, in an approx- 
imate way, the important vacuum-ultraviolet contributions of bound-free and line transi- 
tions. This model was used to obtain solutions pertinent to the case of reentry into the 
earth T s atmosphere. While the results are restricted to small values of the radiation 
cooling parameter, which characterizes the relative importance of radiation and convec- 
tion as energy-transport mechanisms, they cover broad ranges of vehicle velocity, altitude, 
and nose radius, which are of practical interest. 

The characteristic enthalpy variation of the model absorption coefficient was found 
to be nearly independent of altitude and nose radius for fixed vehicle velocity except for 
velocities lower than 10.67 km/sec. Thus it was possible to correlate certain quantities 
by plotting these quantities as functions of the nondimensional adiabatic radiant heat- 
transfer rate for various altitudes and nose radii at fixed vehicle velocity. Among the 
quantities correlated was the cooling factor (the ratio of the stagnation-point radiant heat- 
transfer rate to the adiabatic radiant heat-transfer rate). The cooling-factor correlation 
is particularly useful because it eliminates the need to perform nonadiabatic calculations 
whenever radiant heat-transfer rates are desired. Also correlated was the factor by 
which the convective heat-transfer rate is reduced because of radiation losses in the 
shock layer. Finally, upper-bound estimates were made of the effects of absorption of 
precursor radiation by the free-stream air on the radiant and convective heat-transfer 
rates. 


INTRODUCTION 


At earth entry speeds near, or in excess of, escape speed (about 11 km/sec) blunt 
bodies experience significant heating by thermal radiation from the hot-air shock layer. 

It has been shown that at these flight conditions, neither the influence of energy transport 
by radiation on the thermodynamic and flow properties of the shock-layer gas nor the non- 
gray character of the absorption coefficient of the gas can be ignored if reasonable quan- 
titative estimates of radiant heat-transfer rates are to be obtained. (See, for example, 
refs. 1, 2, and 3.) 

Even though these heating rates are large (from the point of view of heat protection), 
the rate at which energy is transported by radiation is small compared with the rate at 
which energy is transported by the flow for a broad range of conditions of practical inter- 
est. As a result the flow properties are only slightly perturbed from the radiationless 
case. Goulard (ref. 4) took advantage of this fact and developed what amounted to a first- 
order perturbation solution of the temperature distribution in the inviscid stagnation 
region of an optically thin, gray, gas layer. This first-order result corresponds to the 
case of a gas which radiates at a constant rate regardless of its local thermodynamic 
state. 

A second-order perturbation solution was obtained in reference 2. This solution 
includes the effects of temperature -varying thermodynamic and optical properties, which 
are very important because of the strong temperature dependence of these properties. 
However, it was found that the conventional perturbation expansion diverged in the vicinity 
of the body surface, and this fact indicates that the problem is of the singular perturba- 
tion type. The Poincare -Lighthill- Kuo (P-L-K), or perturbation-of-coordinates, proce- 
dure was used to obtain a uniformly valid perturbation expansion. This solution was then 
used to study the effects of various parameters for an idealized gray gas. These results 
(originally presented in ref. 2) are presented herein to provide some physical insight into 
the nature of the influence of the several parameters. 

A nongray absorption-coefficient model based on recent published information on 
the radiative properties of high-temperature air is presented. This model, which includes 
the important features of the vacuum-ultraviolet bound-free continuum and line contribu- 
tions, was used in conjunction with the singular perturbation solution to calculate the flow- 
field properties and stagnation-point radiant heat-transfer rates over a range of free- 
stream velocities (7.92 to 15.24 km/sec), altitudes (30.48 to 73.15 km), and nose radii 
(1 to 1000 cm). A cooling factor, which represents the decrease in radiant heat-transfer 
rate due to radiation cooling of the shock-layer gas, is defined and a correlation in terms 
of the adiabatic (or isothermal slab) radiant heat-transfer rate is presented. A correla- 
tion of the approximate decrease in convective heat-transfer rate due to radiation cooling 
is also presented. 
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The effects of absorption of precursor radiation and of absorption in an ablated 
vapor layer were not included in the calculations. However, if emission is small com- 
pared with absorption in both the free stream and ablated vapor layers (a physically rea- 
sonable assumption because of the relatively low temperatures), the effects of these phe- 
nomena can be determined separately and the results of this paper modified by simple 
correction factors. Upper bounds to the correction factors for the increase in radiant 
and convective heat-transfer rates due to absorption of the precursor radiation are pre- 
sented. No attempt was made to determine the radiation-blocking effect of the ablated 
vapor layer. 


SYMBOLS 


a constant in equation (2) 


B nondimensional blackbody emissive power, T 4 /t s 4 


Bi integral of nondimensional Planck function for ith step 


nondimensional Planck distribution function, Planck distribution functio n 

- T s 4 

7 T s 

^2 

B\ dX 

XI 


Bx 

B(Xl,X 2 ) = 


Cp s specific heat at constant pressure immediately behind shock, cm2/sec2-°K 

E n exponential integral function of order n; E n (y) = \ e - ^t“ n dt 

J 1 

F function defined by equation (A60) 

T-. , , QR/w 

F c cooling factor, — 

^RjWja 

f nondimensional stream function 

f n nth-order coefficient in perturbation expansion of f 

f n * nth-order coefficient in P-L-K expansion of f 

H function defined by equation (A64) 
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h nondime nsional enthalpy, j foth^Py 

— W ^ 

2 W °° 

h e value of h at the edge of the viscous boundary layer 

h n nth-order coefficient in perturbation expansion of h 

h n * nth-order coefficient in P-L-K expansion of h 

l(M>A- 2 >Li) radiant flux for wavelength interval (xi,X2) emitted from one side of infinite 
slab of thickness L, W/cm2 

I x w monochromatic radiant flux absorbed by wall (stagnation point), W/cm3-sr 

kp Bouguer number, p g Kp g A a 

k s thermal conductivity immediately behind shock, W/cm-°K 


k x = k P K X 


L 

N Pe 

q C,w 


slab thickness, cm 

Peclet number (product of Reynolds number and Prandtl number), 

PooWooCp^gAa/kg 

nondimensional convective heat flux at wall (stagnation point) , 
Convective heat flux 



q R f 


nondimensional divergence of radiant-heat-flux vector, 
Divergence of radiant-heat -flux v e ctor 
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co 


q^' function defined by equation (A65) 



nth-order coefficient in the perturbation expansion of q^’ 
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q R,s 


q R, 


w 


nondimensional radiant heat flux leaving front of shock layer in upstream 
direction Radiant heat flux leaving front of shock layer 

\ PooWoo 3 

nondimensional radiant heat flux absorbed by wall (stagnation point), 
Radiant heat flux absorbed by wall 

ip w 3 

9 H oo ■» 


q R,w,a 

R n 


Re 


nondimensional adiabatic radiant heat flux absorbed by wall (stagnation point) 
nose radius in vicinity of stagnation streamline, cm 
shock radius in vicinity of stagnation streamline, cm 
r radial coordinate, cm 

r Q ^,r^ ^ quantities defined by equations (16) and (21), respectively 
coefficient of surface reflectivity 
exponent in equation (46) 


l w 


T 

T c 


W„ 


temperature, °K 

temperature immediately behind normal shock, °K 
radial component of velocity, cm/sec 
free-stream velocity, cm/sec 


W„ 


effective free-stream velocity, Wooj/l + yqj^ )S , cm/sec 


w 


axial component of velocity, cm/sec 


transformed nondimensional Dorodnitsyn coordinate 


axial coordinate, cm 
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0-1 ,0-2 


0 


y 

A 

A. 

A 


e 


C 

V 




^A 

"A,: 


0 l> e 2 

K P 


K 


P,S 


step heights in absorption-coefficient model for line contribution 

relative density of gray background in absorption-coefficient model for line 
contribution 

fraction of precursor radiation absorbed by free stream in path of body 
shock standoff distance, cm 

shock standoff distance for nonradiating (adiabatic) shock layer, cm 

ratio of shock standoff distance for radiating and nonradiating shock layer, 
A/A a 


radiation cooling parameter, 


4oT s 4 kp 


stretched boundary-layer coordinate defined by equation (A61) 
nondime nsional Dorodnitsyn coordinate 
thickness of viscous boundary layer in terms of rj 
shock location in terms of Tj 

nth-order coefficients in perturbation expansion of rj ^ 

nth-order coefficient in P-L-K expansion of rj 

quantities defined by equations (A39) and (A40), respectively 

nondime nsional Planck mean mass absorption coefficient, 

Planck mean mass absorption coefficient 
K P,s 

Planck mean mass absorption coefficient at conditions immediately behind 
normal shock, cm2/g 


6 


nondimensional monochromatic mass absorption coefficient, 
Monochromatic mass absorption coefficient 
*P,s 

X wavelength, A 

Xj,X2 endpoints of wavelength interval, A 

/T Planck mean linear absorption coefficient, cm-1 

jUj linear absorption coefficient for ith step, cm-1 

£ variable of integration 


nondimensional density, P ens ^y 

P s 

standard sea-level density of air, 1.225 x 10"^ g/cm^ 
density immediately behind normal shock, g/cm^ 
free-stream density, g/cm3 


p * 

^oo 


PK 


a 


effective free-stream density, p ( ^°° | , g/cm3 

\w * / 


(M’^2) volume absorption coefficient averaged over wavelength interval 
cm~l 

Stefan-Boltzmann constant, 5.6697 x 10“ ^ W/cm2-°K^ 

normalized optical path length in a gray gas, Optics path length 

k p 

normalized monochromatic optical path length, Qpti ca -1 P at ^ lengt h 




factor for increase of convective heat-transfer rate when absorption of 
precursor radiation by free stream is taken into account 
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<£r factor for increase of radiant heat-transfer rate when absorption of precursor 

radiation by free stream is taken into account 

functions defined by equations (17) and (22), respectively 
<Pl,(j>2 quantities defined by equations (A3 7) and (A38), respectively 
X density ratio across normal shock, P^jpg 

4 / i>4 / 2 quantities defined by equations (A41) and (A42), respectively 

Dots indicate derivatives with respect to hp. 

Primes indicate derivatives with respect to tj. 

ANALYSIS 

Stagnation-Region Flow Model 

Available numerical solutions indicate that a reduction from a three-dimensional to 
a nearly equivalent one-dimensional problem can be carried out in the stagnation region. 
(See, for example, ref. 5.) This simplification can be applied because the flow behind a 
strong bow shock is nearly incompressible in the stagnation region. Also, the various 
thermodynamic properties are nearly independent of the lateral or radial coordinate. 

While the same arguments apply in the stagnation region of a radiating shock layer, it is 
not possible to postulate the existence (even approximately) of a one-dimensional solution 
solely on this basis. Some additional assumption is required regarding the effect of the 
far field on the radiant heat flux and its divergence. This effect is not known a priori as 
it depends on the solution to the entire flow field. Fortunately, the shock layer is thin and 
only a small portion of the radiant energy emitted by gas removed from the stagnation 
region actually passes through the stagnation region. If absorption is small, only a small 
portion of this energy is absorbed in the stagnation region. On the other hand, if absorp- 
tion is large, the beam is greatly attenuated when it reaches the stagnation region, and 
only a small portion of the energy which started the journey is absorbed there. The 
divergence of the radiant flux is influenced only by the amount of energy absorbed and 
emitted. Consequently, the far -field effect on the divergence of the radiant flux is a result 
of that small portion of radiant energy originating in the far field and absorbed in the stag- 
nation region. In the transparent and optically thick limits, this effect of the far field 
vanishes. 

It would appear from the foregoing discussion that a stagnation model for a radiating 
gas can be postulated as long as the assumptions concerning the far field are not grossly 
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unrealistic. In what follows, a particular stagnation 
model will be formulated and an estimate of the inac- 
curacy resulting from the assumption concerning the 
far field will be obtained. 

A schematic diagram of the flow in the stagna- 
tion region of a blunt body is shown in figure 1. At 
very high speeds, the ratio of the shock standoff dis- 
tance A to the shock radius R s is very much 
smaller than 1 (a typical value is 0.05). Under these 
conditions, the geometry of the stagnation region 
closely resembles a plane parallel gas slab. In addi- 
tion, the enthalpy in the shock layer varies slowly with 
respect to r/A so that the stagnation region may be 
approximately represented by a gas slab in which the 
thermodynamics as well as the geometry is one 
dimensional. 


Figure 1.- Schematic digram of flow As & resu ^ 0 f these considerations, the model 

described herein has been chosen to represent the 
flow in the stagnation region of a blunt object. The model consists of an axially symmet- 
ric flow impinging upon a flat plate of infinite extent oriented so that the normal to the 
plate is parallel with the stream direction. At a plane which is parallel with the plate and 
at a distance A in front of it, the gas is suddenly raised to a total specific enthalpy of 
1 2 

2 Woo . The geometry of the flow model is illustrated in figure 2. 


The boundary conditions on the radiant- energy flux are specified as follows: 


(1) The boundary at z = A (which corresponds to a bow shock) is transparent. 


(2) There is no radiant-energy transfer from the free 
stream to the shock layer. 

(3) The boundary at z = 0 (which corresponds to the 
body surface) is cold and reflects diffusely and indepen- 
dently of wavelength a fraction r w of the incident 
radiation. 

The statement that the body surface is cold means 
that emission from the body surface has a negligible influ- 
ence on the gas in-the shock layer. When the hot shock 
layer (temperatures in excess of 10 000° K) is optically 
thin, emission from the relatively cool body surface (tem- 
peratures less than 4000° K) may be comparable to 



\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 
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Figure 2.- Geometry of the plane- 
parallel-slab model. 
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emission from the shock-layer gas. However, because the shock layer is optically thin, 
very little of the radiant energy emitted by the body surface will be absorbed by the shock- 
layer gas. On the other hand, when the shock layer is optically thick and absorption is 
important, the shock-layer gas emission will approach the blackbody value corresponding 
to the high temperature of the shock layer. Since blackbody radiation is proportional to 
the fourth power of temperature, the energy emitted from an optically thick gas layer will 
greatly exceed that emitted from the body surface. Thus, whenever the body surface tem- 
perature is small compared with the shock-layer temperature, the influence of emission 
from the body surface on the shock-layer gas is unimportant. 

Governing Equations 

The governing nondimensional equations for the proposed inviscid stagnation flow 
model are 


f(77)h’(rj) + q R '(r?) = 0 (1) 

2ffa)f"(Tj) - [f'(7?)] 2 + ^ = 0 (2) 

The boundary conditions are 

f(0) = 0 (3) 

f^A) = 1 (4) 

< 5 > 

\ I xRs l/2x(l - X) 

h (t?a) = 1 (6) 


The quantity p g A a ?7 is the Dorodnitsyn coordinate, 2p 0O W co f(?7) the stream function, 

■i the total enthalpy, PooW^q ’(77) the divergence of the radiant heat-flux 

vector, A a the adiabatic (or radiationless) shock standoff distance, R s the shock 
radius, x = P<x>/Ps density ratio across the shock, and p s p the density. The sub- 
script s has been used to denote dimensional quantities immediately behind the normal 
shock. The quantity p g A a ? 7 ^ defines the location of the shock in terms of the 
Dorodnitsyn coordinate. 


10 



This system of equations is similar to the system obtained by Howe and Viegas 
(ref. 6) when viscosity, thermal conductivity, and curvature effects are neglected. A 
complete derivation of the system of equations (1) to (6) is presented in reference 2. 
Implicit in this derivation is the strong- shock assumption, which leads to the neglect of 
the normal pressure gradient. The shock standoff distance (and, hence, the constant a) 
is related to the lateral pressure gradient and the density ratio through the Rankine- 
Hugoniot conditions across a spherical shock. 

The nondimensional form of the divergence of the radiant flux vector was also 
derived in reference 2: 


q R '(v) = -2e/cp (77)13(77) + ekpj^ A K x (£)B x (£)Ei^kp 7^(77) - r^) 

+ 2r w E 2 ^kpr x (r ? ))yJ A K x («)B x (|)E 2 (kpT x (*))d* 


m 


(7) 


The quantity /<p g Kp is the Planck mean mass absorption coefficient, Kp g K x the 

monochromatic mass absorption coefficient, oT s ^B^ 7 r the blackbody emissive power, 
oT s ^B^jn the Planck function, and E n the exponential integral functions of order n. 
The optical path length is 

rV 

k pTx^) = k P \ P K \ d7 7 ( 8 ) 

It is the presence of the second and third terms on the right-hand side of equa- 
tion (7) which so greatly complicates the radiation problem. These terms are integral 
expressions. In addition, their presence makes it impossible to define a wavelength- 
averaged absorption coefficient by which the wavelength dependence might be eliminated. 
The importance of these terms is indicated by the magnitude of the Bouquer number, 
k P = Ps K p s ^a> ratio of the shock standoff distance for a nonradiating shock 

layer to the photon mean free path evaluated at conditions immediately behind the shock. 

The radiation cooling parameter, e = (^oT s ^ kp, admits of several physical 
interpretations which are useful in the understanding of the radiating shock layer. Of 
these, one of the most useful is the following: 

I Rate of emission from element of yTime required by element of volume to traverse! 
I volume of gas emerging from shock) (distance A a at rate of emergence from shock J 

2(Energy of element of volume upon emergence from shock) 
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This parameter modifies the entire radiation term, and thus it is a measure of the effi- 
ciency of radiation compared with convection as an energy transport mechanism within 
the shock layer. In addition, the surface reflectivity r w and the wavelength dependence 
of the absorption coefficient influence the character of the radiation terms and will be 
considered as parameters in this study. 

A relationship between the thermodynamic variables h and p is needed to com- 
plete the set of equations. The simple correlation formula (obtained in ref. 2 ) 



was used in this study. This formula is generally accurate to about 5 percent. 

Conventional Perturbation Procedure 

When the radiation cooling parameter e is very small, the radiation terms in the 
energy equation (eq. (1)) become of only secondary significance throughout most of the 
domain of the problem. 1 Neglecting the term qp ^(77) reduces the problem to one in 
which radiation plays no part. If, as expected when e is small, the presence of radia- 
tion only slightly influences the solution, qp'(ry) can be evaluated with reasonable accu- 
racy by using the radiationless solution for h so that equation (1) becomes purely dif- 
ferential. Then, when the small-perturbation procedure (which is carried out approxi- 
mately in this manner) is applied to this problem, the integrodifferential system is 
simplified to a purely differential system. In addition, as a result of the nature of the 
radiationless solution for the enthalpy distribution, the energy and momentum equations 
become uncoupled and can be solved independently. Hence, analytic solutions to the flow 
in the inviscid region of the shock layer can be obtained to any order of approximation. 
Details of the derivation of these solutions are presented in appendix A. 

The zero-order, or radiationless, solution is simply 


h 0 (?7) = 1 


( 9 ) 


fo fa) = (1 - a)?? 2 + ary 


( 10 ) 


where 


2\/2x(l - X) 

1 + \/2x(l - x) 


( 11 ) 


■^With the obvious exception of the region 77 ~ 0 where f (77) ~ 0. 
presented by this exception are discussed subsequently. 


The difficulties 
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The first-order solution, which represents the effect of radiation when the optical and 
thermodynamic properties of the gas are assumed to be independent of the enthalpy, is 


fl 

h l(?l) = \ 2 o 

J ri (1 - a + a£ 


h(v) = - it - a)„ - ajj'J - lK> d{ ♦ r? f 2(1 - - a ^ 

[(1 - a)ij + aj ij 2 |(l - a)£ + a] 


( 12 ) 


$lU)dO (13) 


Here £, is a dummy variable of integration. The quantity q-^ q'(v) is given by the 
formula 


aoo r 

%R,o'W = -\ 0 K X B X E 2( k xd - -n)) + (l - r 0)X )E 2 (k x 7 7 ) 


dX 


(14) 


The notation has been simplified somewhat in this expression by omitting the argument 
hq in the terms k-^ and B\ and by introducing the quantities 


Also 


k X = k P*X 

r 0,X = r w[i - 2E 3( k x) 


$l(77) = a 2 h!(rj) 


(15) 

(16) 

(17) 


The second-order solution takes into account the change in gas properties with 
changes in enthalpy. This solution is 


h2(?7) = 7 A> iq R>0 , (l) + J 


q R,l t ^ f l^)^R,O t ^) 


fo(€) 


f 0 2 (^) 




(18) 


f2^) = "(I - a ) 2? 7A,l 27?2 " ^j[ 2 (l _ + 2 ^ 


f 

a) | + a]' 


d« 


+ r? 2 C 


2(1 - a)| + a 


n l 2 [(l - a ) 4 + a]“ 




(19) 
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where 


<lR,l'fa) = h l ( ^I 0 |^A B x[ E 2(kx(l - V)) + (l - r 0>X ) E 2 (k x ?7) 

+ C h l ( ^ )E l( k ^ 77 - *l) d * + k 


~ 2 «A B X> dX 


1 a b a 


L (k x d 


(l-r 0 ,A .) E lM)C h l(« d « 


+ r l,A E 2( k A^)J> d ^ 


( 20 ) 


Here the argument hQ is omitted in the terms k x and B x and the quantity rj x is 
defined by the expression 


r l,A = 2r W 


Also 


kxBaJ 1 hiU)E 2 (k x £)d* + A x B x E 2 (k x ) hi($)d* + r? A>1 ^B x E 2 (k x 


2 

*2fo) = -fi(T7)fi ,r (T?) + - | a 2 h 2 (r?) 


The quantities ? 7 A q, r? A j, and rj A 2 are given by the formulas 


( 21 ) 


(22) 


^,0 


= 1 


*?A 


-If 1 

2 J n 


*lU) 


0 [(1 - a) | + aj 


d| 


9 1 pi 

^,2 = ^ - ^A,! + 2 J 0 




[(1 - a) | + a] 


d£ 


(23) 


(24) 


(25) 


It can be seen upon inspection of equation (20) that a large rate of change of the 
absorption coefficient with enthalpy will lead to large values of q^ Thus, it is 

clear that for conditions of interest, for which the absorption coefficient does vary strongly 
with enthalpy, the second-order solutions can become more important to the overall 
solution than their order in 6 might at first indicate. 
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Perturbation-of-Coordinates Procedure 

As can be seen from an inspection of the expressions (12) and (18), the first-order 
solution for the enthalpy distribution has a logarithmic singularity at the point r\ = 0 and 
the second-order solution has a singularity of greater strength at this point. As a con- 
sequence, the assumed expansion diverges as the origin is approached and the small- 
perturbation solution is not uniformly valid throughout the domain of the problem. This 
divergence can lead to serious errors in the calculation of the radiant heat flux to the wall 
because those regions close to the wall, in which the largest errors occur, are given the 
most weight in the calculation. This is particularly true for shock layers that are not 
optically thin. Additional difficulties are encountered when attempting to specify the 
proper outer boundary conditions for the viscous boundary -layer equations. In classical 
boundary -layer theory, the outer boundary conditions are obtained from the values of the 
outer (or inviscid) solution at the wall (77 = 0 in this problem). Because of the divergence 
of the outer solution, no finite value exists at rj = 0. 

In this section the method of perturbation of coordinates^ (ref. 7) is used to obtain 
a solution which is uniformly valid over the domain of the problem. Details of the appli- 
cation of this method to the problem of this paper are presented in appendix A. This 
method uses a coordinate transformation in the form of a perturbation expansion of the 
coordinate to remove the singularity (which caused the divergence of the conventional 
solution) from 77 = 0 to a small negative value of 77 which lies outside the domain of 
the problem. The perturbation-of-coordinate expansions are 

7 } = x + e?7^*(x) + . . . (26) 

h(r],e) = h 0 *(x) + ehi*(x) + . . . (27) 

f(T 7 ,e) = f 0 *(x) + efi*(x) + . . . (28) 

where x is the coordinate in the transformed plane, and the asterisk has been used to 
differentiate between the coefficients in the perturbation-of-coordinate expansion and the 
coefficients in the conventional expansion. Pritulo (ref. 8 ) has derived a general relation 
between the coefficients in the two types of expansion. Adapted to this problem, the rela- 
tionships became 


h 0 *(x) = h 0 (x) 

(29) 

hi*(x) = hj(x) 

(30) 


^Variously called the Poincare-Lighthill-Kuo method, the P-L-K method, the P-L 
method, LighthilPs technique, and the method of strained coordinates. 
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f 0 *(x) = f 0 (x) 

(31) 

fl*(x) = fl(x) + 77j*(x)fo’(x) 

(32) 

*t \ h 2 (x) 

?7l*(x) — 

hi (x) 

(33) 


The second-order term h2(x) introduces the effects of variable thermodynamic and 
optical properties, so it is apparent that these effects are contained in the first-order 
perturbation-of-coordinate solution. 

A comparison of the perturbation-of-coordinate and conventional perturbation solu- 
tions for the enthalpy distribution for a constant-density, transparent shock layer is pre- 
sented in figure 3(a). The divergent character of the conventional solutions is apparent. 



(a) Enthalpy distribution, e = 0.1; kp = 0; (b) Percent error in radiant flux, 

p = Constant; KpB = h 6 . 

Figure 3.- Comparison of the P-L-K and conventional perturbation solutions. 
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Also shown in the figure is the exact analytic solution, which can be found in this simple 
case. The formula for this exact solution is 


h(r?) = 


/I + 4e(d - l)ln 

1 + 

l 

_ 2 (’ J /”a) Jj 


1 

d-1 


(34) 


where d (the exponent in the correlation formula KpB = h^) was taken to be 6 and the 
constant a (which appears in the momentum equation, eq. (2)) was taken to be 0.5. The 
good agreement between the perturbation-of-coordinate solution and the exact solution 
indicates that the accuracy of the perturbation-of-coordinate solution is probably second- 
order in the radiation cooling parameter e throughout the domain, except perhaps in the 
immediate neighborhood of the wall. It is clear that quantities such as the radiant heat 
flux at the wall, which depend upon an integration over the enthalpy distribution, will be 
considerably more accurate if the perturbation-of-coordinate solution rather than the con- 
ventional perturbation solution is used. 

The percent error in radiant heat flux is presented as a function of e in figure 3(b). 
The error is less than 12 percent for values of e no larger than 0.3. 

It should be noticed that the perturbation-of-coordinate solution does not lead to 
zero enthalpy at the wall as the exact transparent solution does. The reason for this 
disparity can be found in the fact that the coordinate stretching displaces the boundary 
with regard to both the energy and momentum equations but not by a uniform amount. 

Thus a physical interpretation of the first-order perturbation-of-coordinate solution is 
that the normal velocity of the flow at the boundary for the energy equation is not quite 
zero, and a particle approaching this boundary will reach it in a finite time before losing 
all its energy by radiation. 

Since the expected error in the Dorodnitsyn coordinate 7] in terms of the stretched 
coordinate x is of order e2 and since the gradients in h^(x) are very large in the 
vicinity of the wall, the difference between the perturbation-of-coordinate and exact solu- 
tions at the wall lies within expected limits. Convergence to the correct solution should 
be attained with the addition of higher order terms to the expansions of h and rj. 


RESULTS AND DISCUSSION 


Results of certain of the calculations in this report showed that the enthalpy is a 
double- or even triple-valued function of the Dorodnitsyn coordinate r) in the vicinity of 
the shock for large values of the Bouguer number. (See fig. 4 for typical example.) An 
examination of the governing equations failed to show the presence of any singularities 
that might adversely influence the solution in this region when kp is large and € is 
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small. On the other hand, the results of calculations with varying mesh size seemed to 
rule out the possibility that this physically unrealistic behavior can be attributed solely 
to numerical inaccuracies. Consequently, it is suspected that the difficulty results from 
the large first and second derivatives of h{rj) just behind the shock, which can lead to 


1.0 r 

.8 “ 

.6 - 
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Figure 4.- Example of triple-valued solution, e = 0.1; Kp = 10; Kp =4; r w = 0. 

values of T7j*(x) very much larger than 1. This behavior is to be expected because of 
the boundary-layer type of behavior of the enthalpy just behind the shock for an optically 
thick shock layer (see ref. 2). It is suspected that inclusion of higher order terms in 
expansions would either eliminate the problem or increase the value of kp at which it 
first appears. For truncation after the second-order terms the condition for validity of 
the solution appears to be ekp < 1. 


Gray -Gas Results 

Shock-layer enthalpy distributions for a gray gas with differing values of the radia- 
tion cooling parameter e, the Bouguer number kp, the variation with enthalpy of the 
Planck mean mass absorption coefficient k-p, and the reflectivity of the body surface r w 
are presented in figures 5 to 7. While the gray-gas assumption may not be realistic for 
most gases of interest, its use is felt to be justified in the study of these parameters for 
two reasons. First, the highly complex and varied spectral structure of absorption coef- 
ficients makes a general parametric study of nongray gases impractical. Second, experi- 
ence with nongray calculations indicates that the qualitative dependence of the gray results 
on the various parameters will carry over to most nongray cases. 
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The value of x> the density ratio across the normal shock, was fixed at 0.06 for 
the calculations reported on in this section. This choice is justified because x varies 
but little with altitude and velocity and the effects of this variation on the stagnation solu- 
tions are slight. The value x = 0*06 is typical for hypervelocity flight in the atmo- 
sphere of the earth. 



(Ct £ = 0.10. 

Figure 5.- Concluded. 

The decrease in enthalpy level with increasing e is illustrated in figures 5(a) to 
5(c). These results indicate that the loss of energy from the shock layer by radiation 
(i.e., radiation cooling) can produce a noticeable drop in enthalpy for values of e as 
small as 0.01. The dependence of the enthalpy distribution on the Bouguer number (hence, 
optical thickness) is also shown in these figures. As expected, an increase in the Bouguer 
number (or optical thickness) inhibits shock-layer cooling and leads to higher values of 
enthalpy near the wall. 

The variation of the enthalpy distribution with Kp (the enthalpy variation of the 
Planck mean mass absorption coefficient) for several values of the Bouguer number kp 

is shown in figures 6(a) to 6(c). These effects are most noticeable for optically thin shock 
layers (kp « l.o) and tend to vanish as the optical thickness increases. In a transparent 
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layer, the rate of emission of radiant energy is proportional to the Planck mean mass 
absorption coefficient Kp. Thus, gases with small values of k p (which mean larger 
values of Kp when the nondimensional enthalpy is less than 1) will be cooled more than 
gases with large values of kp. As the optical thickness increases, a smaller value of 
kp still implies greater emission rates but it also means greater absorption and more 



(c) kp = 1.0 
Figure 6.- Concluded. 

radiant energy available for absorption. The process of absorption tends to counteract 
the differences in emission rates due to differences in kp. Finally, when radiation equi- 
librium is reached (this state is achieved in the interior of optically thick regions), the 
energy of the particle is independent of its optical properties. Of course, in those regions 
optically close to the shock and the wall the amount of radiant energy available for absorp- 
tion is not so great as in the interior of the shock layer, and particles in these regions 
cannot approach the state of radiation equilibrium (except in a region optically close to a 
perfectly reflecting surface). Thus, the enthalpy distribution remains dependent on the 
value of kp near the shock and the wall. This dependence of kp is suppressed near 
the shock, where h is almost 1, because the values of Kp are nearly the same despite 
the differences in Kp. 
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The effect of surface reflectivity r w on the shock- layer enthalpy distribution is 
shown in figures 7(a) to 7(c). If the shock-layer gas is transparent (i.e., the gas does 
not absorb) , surface reflectivity has no effect on the enthalpy distribution because all 
photons emitted by the layer escape. Whether or not a photon is absorbed or reflected 
by the wall is of no consequence. As the optical thickness of the layer is increased, the 



(a) k P = 0.1. 

Figure 7.- Effect of surface reflectivity on the shock-layer enthalpy 
distribution. £ = 0.1; Kp = 4.0. 

chance of capture of a photon by absorption in the layer is increased. If the surface 
reflectivity is increased also, the probability of capture is increased still further because 
many photons which might otherwise have escaped into the wall are reflected back into the 
layer and are once again subject to capture there. Consequently, the enthalpy level is 
higher near a reflecting wall than it would be near a nonreflecting wall. 

It can be concluded from the foregoing discussion that use of a reflecting surface 
will not reduce the radiant heat-transfer rate from the gas to the wall by the factor 
1 - r w (unless, of. course, the gas is transparent) but will reduce it by some smaller 
fraction. The reason is that the radiant heat flux incident on the wall is larger when the 
wall is reflecting because of the higher enthalpy level. In addition, the rate of heat trans- 
ferred to the wall by conduction will be greater, also because of the higher enthalpy level. 
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Of course, increasing the surface reflectivity always decreases the total heat-transfer 
rate to the wall because the higher enthalpy level must lead to an increased loss of energy 
by radiation through the shock in the upstream direction and by convection in a lateral 
direction away from the stagnation point. 

Figures 8 to 10 show the effects of variations in the parameters e, kp, Kp, and 

r w on w , the rate of radiant heat transfer to the wall normalized by i Poo , the 

energy influx to the shock layer. The rate of radiant heat transfer to the stagnation point 
was calculated from the formula 

^R,w = e ( x “ K p( ? l) B ( T l) E 2 ( K p T (v)J d V (35) 

where the optical thickness kpT is given by 


rV 

k p T(?]) = k PJ Q K p (r])dr] 


(36) 


The dashed curves in figure 8 indicate the "no-decay limits" for various values of 
the Bouguer number. These limiting curves were computed by assuming the shock layer 

to be isenthalpic so that Kp{r}) = B(r;) 
= 1. Thus, the no-decay limit curves 
are given by the formula 


P-L-K solution 
Isenthalpic 



Figure 8.- 


Effect of £ and kp on the rate of radiant heat transfer 
to the stagnation point. Kp = 4.0 ; r w = 0. 


q R, 


w 


= e (! - r w) 


1 - 2E, 


;( k p) 


2k T 


(37) 


where E3^kp) is the exponential 
integral function of third order. This 
no-decay approximation is often used 
to predict the rate of radiant heat 
transfer when radiation effects are 
small. Use of this approximation 
always gives an upper bound to the 
true value of qp w * A study of fig- 
ure 8 indicates that the no-decay limit 
curve is least accurate in predicting 
the rate of radiant heat transfer in the 
transparent case kp = 0. This result 
is expected because the enthalpy 
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Figure 9.- Effect of the enthalpy variation of the absorption 
coefficient on the rate of radiant heat transfer to the 
stagnation point, e = 0.1; r w = 0. 


in which Kp is least, 
by differences in the value 
increases. 


distribution for the transparent case is the 
most perturbed from an isenthalpic state. 
Results presented in this figure also indi- 
cate the importance of absorption (as char- 
acterized by the Bouguer number kp) in 

reducing the rate of radiant heat transfer 
from the shock layer to the wall. 

The results presented in figure 9 
indicate that the differences in /ip are 

most important when the optical thickness 
of the shock layer is small. Here the radi- 
ant heat transfer to the wall is greatest for 
the smallest value of /ip. This, of course, 


supplements the observation (from fig. 6(a)) 
_ that radiation cooling is greatest for gases 
The differences in radiant heat transfer to the wall brought about 
of /ip tend to vanish as the optical thickness of the layer 


The reduction in radiant heat transfer to the wall due to surface reflectivity is 
shown in figure 10. When the shock layer is transparent, the rate of radiant heat trans- 
ferred q R is in the ratio 1 - r w . However, as the optical thickness of the shock 

layer increases, the ratio becomes somewhat greater than 1 - r w as predicted earlier 
in this section. 


The effect of the parameters 
e, Kp, and r w on the shock 

standoff distance is shown in fig- 
ures 11 and 12. The quantity A 
is the ratio of the shock standoff 
distance in a radiating shock 
layer to that in a nonradiating 
(or adiabatic) shock layer at the 
same flight conditions. It was 
computed from the formula 



A = 



(38) 


Figure 10.- Effect of surface reflectivity on the rate of radiant heat 
transfer to the stagnation point, e = 0.1; Kp = 4.0. 
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Figure 11.- Effect of e and kp on the shock standoff distance, 
kp = 4.0; r w = 0. 
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Figure 12.- Effect of surface reflectivity on the shock standoff distance. 

£ = 0.1; kp = 4.0. 

The results shown in figures 11 and 12 indicate, as expected, that a decrease in 
enthalpy level (with the consequent increase in density level) in a shock layer leads to a 
reduction in shock standoff distance. 

Nongray Absorption-Coefficient Model 

The absorption coefficient of high-temperature air and other gases and gas mixtures 
depends strongly on wavelength. Consequently, the gray-gas assumption (i.e., that the 
optical properties of the gas are independent of wavelength) is poor indeed, and has been 
resorted to so frequently in the literature only because of the resulting simplicity. The 
small-perturbation solution derived in this paper overcomes these difficulties by reducing 
the absorption integrals in the divergence of the radiant flux vector to a form amenable to 
direct evaluation. Thus, it is necessary only to perform an integration over a known, 
albeit complicated, function of wavelength. In view of the current uncertainties with 
regard to spectral distributions of gaseous absorption coefficients, a simplified step- 
function model was used for the absorption coefficient for air. 
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Conditions of interest in this paper are the temperature range of 8000° to 16 000° K 
and the density range of 10“^ to 10“* times standard sea-level density. The upper limit 
on the temperature was chosen so that the arbitrarily set condition e < 0.3 for the valid- 
ity of the small-perturbation solution would not be violated. The range of altitudes and 
free-stream velocities which produce these conditions behind a normal shock are shown 
in figure 13. At these conditions radiation can be an important mechanism for energy 
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Figure 13.- Range of conditions investigated, shown by the shaded area. 


transport. The gas in the stagnation region of the shock layer generally may be consid- 
ered to be in chemical equilibrium (this, of course, also depends on body size), and the 
small-perturbation solution of this paper is generally valid (this also depends on body 
size). Several typical entry trajectories have been drawn on the altitude -velocity chart of 
figure 13. It is apparent that the range of conditions to be studied herein is representative 
of the environment encountered by a large class of entry vehicles. 

Also shown in this figure is the range of altitudes and velocities covered by the cal- 
culations to be discussed. Limitations on this range other than those imposed by the 
limits in the temperature and density are discussed subsequently. 

The spectral distribution of the absorption coefficient of air is characterized by 
large values in the vacuum ultraviolet (wavelengths less than about 1200 A) due to bound- 
free transitions from the low-lying energy states of the nitrogen and oxygen atoms. In 
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addition, large values occur over small wavelength increments because of "line" or 
bound-bound transitions between excited states of the nitrogen and oxygen atoms. The 
strongest of these lines appear in the wavelength region between 900 and 1800 A. In other 
regions of the spectrum where the major contributions are due to bound-free transitions 
from excited states of the nitrogen and oxygen atoms, weak atomic lines, and free-free 
transitions produced by collisions between electrons and the ions of nitrogen and oxygen, 
the value of the absorption coefficient is not so large. . 



Figure 14.- Step-function absorption-coefficient model. 
T = 14 000° K; P S P/P 0 = 10' 2 . 


The step-function model used in this paper is composed of four high steps in the 
vacuum ultraviolet which correspond to the bound-free transitions from the three lowest 
energy levels of the nitrogen atom and the lowest level of the oxygen atom; a set of high, 
narrow, equally spaced steps of uniform height superimposed on a uniform background 
which corresponds to the bound-bound transitions in the ultraviolet; and a single low step 
which covers the rest of the spectrum. The values for the heights of the four vacuum- 
ultraviolet steps were determined by the method described in appendix B from the optical 
properties presented by Hahne in reference 9. The model for the ultraviolet bound-bound 
transitions was used because of the lack of detailed spectral information. Allen (ref. 10) 
presents the data on line contributions in the form of spectrally integrated radiation fluxes 
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emerging from one side of an infinite, isothermal slab of finite thickness (specifically, 

0.1, 1, 5, 10, 15, and 20 cm) for all lines with wavelength less than 2000 A. 

The heights and widths of the tall narrow steps and the height of the background for 
the model absorption coefficient were first obtained by comparing the analytic expression 
for the emergent spectrally integrated radiant fluxes from gas slabs of three different 
thicknesses determined with this model with Allen's results for corresponding slab thick- 
nesses (details are presented in appendix B). However, the resulting values displayed 
some anomalies which were believed to be associated with (1) the inability to read accurate 
values of flux from the curves presented by Allen and (2) the neglect by Allen of the effects 
of line overlapping on self-absorption. Consequently, a program was developed to com- 
pute the flux emergent from isothermal slabs of arbitrary thickness due only to lines in 
the ultraviolet. The parameters used to compute the absorption coefficient for the lines 
were those used by Allen (i.e., those listed by Griem in ref. 11). The expression used to 
compute the line half-widths (see appendix B) was not that used by Allen but was similar 
and gave approximately the same values. The results of this program were then used to 
evaluate the step heights and widths for the model absorption coefficient of the lines in 



Figure 15.- Variation of the radiation cooling parameter with altitude and velocity. 
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TABLE I.- PARAMETERS FOR STEP-FUNCTION MODEL OF ABSORPTION COEFFICIENT 


(a) Linear absorption coefficients for individual steps 


[Final numbers listed with values of pi represent powers of 10; 
for example, 5.001-1 signifies 5.001 X 10”*] 


T, °K 

p aP/P 0 

Ml, 

cm“l 

M 2 , 

cm~l 

M3> 

cm“l 

M4, 

cm~l 

M 5 , 

cm - * 

V 6 > 
cm'l 

M 7> 
cm“l 

M 8 > 

cm"* 

M 9 > 

cm”l 

M, 

cm - * 

M 

8 000 

10 

-3 

5.001-1 

8.629-2 

2.786-2 

8.360-1 

3.475-3 

8.115-1 

3.404-5 

8.082-1 

4.285-6 

7.310-6 

10 000 



4.606-1 

1.067-1 

5.012-2 

2.464+0 

8.511-3 

2.422+0 

2.443-4 

2.414+0 

4.571-5 

1.456-4 

12 000 



3.319-1 

1.004-1 

5.433-2 

1.601+0 

1.132-2 

1.588+0 

5.470-4 

1.547+0 

1.690-4 

6.690-4 

14 000 



1.577-1 

5.902-2 

3.350-2 

9.039-1 

8.091-3 

8.784-1 

6.152-4 

8.710-1 

2.786-4 

1.179-3 

16 000 



5.297-2 

2.291-2 

1.355-2 

4.229-1 

3.681-3 

4.100-1 

3.945-4 

4.097-1 

2.564-4 

1.039-3 

8 000 

10 

-2 

4.678+0 

8.418-1 

2.583-1 

3.165+0 

3.214-2 

2.939+0 

2.570-4 

2.907+0 

8.318-5 

1.152-4 

10 000 



4.880+0 

1.123+0 

5.304-1 

3.029+0 

8.913-2 

2.588+0 

1.242-3 

2.500+0 

4.943-4 

1.322-3 

12 000 



4.389+0 

1.294+0 

7.241-1 

2.853+0 

1.475-1 i 

2.276+0 

3.236-3 

2.132+0 

2.415-3 

8.130-3 

14 000 



3.370+0 

1 . 212+0 

7.295-1 

2.502+0 

1.702-1 ! 

1.943+0 

5.483-3 

1.778+0 

7.194-3 

2.467-2 

16 000 



2.107+0 

8.833-1 

5.480-1 

1.923+0 

1.417-1 

1.516+0 

6.546-3 

1.381+0 

1 . 002-2 

3.958-2 

8 000 

10 

-1 

3.296+1 

7.394+0 

1.717+0 

4.628+0 

2.123-1 

3.124+0 

1.148-3 

2.912+0 

2.254-3 

2.426-3 

10 000 



4.764+1 

1 . 112+1 

5.153+0 

7.941+0 

8.553-1 

3.648+0 

5.888-3 

2.793+0 

6.887-3 

1.350-2 

12 000 



4.775+1 

1.393+1 

7.881+0 

1.050+1 

1.590+0 

4.206+0 

1.507-2 

2.631+0 

2.793-2 

8.247-2 

14 000 



4.375+1 

1.545+1 

9.548+0 

1.189+1 

2,192+0 

4.529+0 

2.594-2 

2.363+0 

8.531-2 

2.939-1 

16 000 



3.715+1 

1.517+1 

9.820+0 

1.189+1 

2.477+0 

4.543+0 

3.412-2 

2 . 100+0 

1.745-1 

6.629-1 


= V/iiBi 

J l ’ ~ l 

(b) Integrals of the nondimensional Planck function* across individual steps 


[Final numbers listed with values of B^ represent powers of 10; 
for example, 1.140-6 signifies 1.140 x 10" 6 ] 


T, °K 

P S P 

/Po 

Bl 

b 2 

B 3 

b 4 

b 5 

B 6 

b 7 

b 8 

b 9 

8 000 

10 

-3 

1.140-6 

2.554-6 

1.841-5 

2.838-9 

7.014-5 

1.079-8 

9.506-3 

1.466-6 

9.901-1 

10 000 



4.123-5 

6.077-5 

3.006-4 

1.429-7 

7.943-4 

3.776-7 

3.811-2 

1.816-5 

9.605-1 

12 000 



4.133-4 

4.410-4 

1.702-3 

1.828-6 

3.524-3 

3.784-6 

8.730-2 

9.397-5 

9.064-1 

14 000 



2.018-3 

1.654-3 

5.358-3 

7.031-6 

9.290-3 

1.219-5 

1.471-1 

1.932-4 

8.343-1 

16 000 



6.342-3 

4.153-3 

1.180-2 

1.315-5 

1.795-2 

1.995-5 

2.055-1 

2.286-4 

7.540-1 

8 000 

10 

-2 

1.140-6 

2.554-6 

1.841-5 

1.047-8 

7.015-5 

3.990-8 

9.506-3 

5.408-6 

9.901-1 

10 000 



4.123-5 

6.077-5 

2.992-4 

9.205-7 

7.925-4 

2.438-6 

3.802-2 

1.170-4 

9.605-1 

12 000 



4.133-4 

4.410-4 

1.690-3 

1.309-5 

3.499-3 

2.710-5 

8.690-2 

6.730-4 

9.064-1 

14 000 



2.018-3 

1.654-3 

5.297-3 

6.668-5 

9.183-3 

1.156-4 

1.455-1 

1.828-3 

8.343-1 

16 000 



6.342-3 

4.153-3 

1.167-2 

1.687-4 

1.770-2 . 

2.564-4 

2.028-1 

2.931-3 

7.540-1 

8 000 

10 

-1 

1.140-6 

2.554-6 

1.837-5 

2.004-7 

6.998-5 , 

5.260-8 

9.484-3 

2.716-5 

9.901-1 

10 000 



4.123-5 

6.077-5 

2.958-4 

1.279-5 

7.816-4 1 

4.831-6 

3.750-2 

6.138-4 

9.605-1 

12 000 



4.133-4 

4.410-4 

1.629-3 

1.517-4 

3.373-3 ! 

7.328-5 

8.375-2 

3.767-3 

9.064-1 

14 000 



2.018-3 

1.654-3 

4.920-3 

7.638-4 

8.551-3 ; 

4.395-4 

1.352-1 

1.208-2 

8.343-1 

16 000 

1 


6.342-3 

4.153-3 j 

1.038-2 

2.198-3 

1.578-2 ! 

1.449-3 

1.806-1 

2.518-2 

7.540-1 


* 

r 852 


r 911 


r 1020 

Bi = 

, B X <& 

b 2 = 

B\ dX 

b 3 = 

P\ B X dX 
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Figure 16.- Variation of the Bouguer number with altitude and velocity. 

the ultraviolet. This model should give a reasonable estimate of the effects of absorp- 
tion, provided that the values of slab thickness chosen do not differ greatly from the 
shock standoff distance and that the temperature distribution in the shock layer does not 
differ much from isothermal. 

The maximum monochromatic absorption coefficient at wavelengths greater than 
1800 A is small enough to insure that the maximum monochromatic optical thickness of a 
shock layer will be very much less than 1 for the temperatures, densities, and shock 
standoff distances of interest. Consequently, the contributions of the various radiators 
can be legitimately accounted for with a Planck type of average absorption coefficient for 
the wavelength interval of 1800 to 100 000 A. 

An example of the model absorption coefficient is shown in figure 14. A compila- 
tion of all the model absorption-coefficient parameters is presented in table I. 

The absorption-coefficient model was used to compute ej R^, the radiation cooling 
parameter per centimeter of nose radius, and kp/R^f, the Bouguer number per centi- 
meter of nose radius, for a range of altitudes and velocities. The results are presented 
in figures 15 and 16. 
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Nongray-Gas Results 


The nongray absorption-coefficient model developed herein was used in conjunction 
with the small-perturbation method to obtain stagnation-point solutions over the range of 
conditions indicated in figure 13. This range of conditions was limited by the tempera- 
ture and density range of the absorption-coefficient model. It was also restricted by 
excluding the high-altitude, low-velocity regime because of large values of the enthalpy 
variation of the Planck mean absorption coefficient at normal-shock conditions, which can 
lead to inaccuracy in the small-perturbation solution. (Fortunately, this is a range in 
which radiation is not particularly important, anyway.) Finally, it was restricted to 
velocities less than or equal to 15.24 km/sec (50 000 ft/sec) because the radiation cooling 
parameter would exceed 0.3 (the arbitrary upper limit for validity of the present method) 
at higher velocities, except for rather small nose radii for which the assumptions of 
thermodynamic and chemical equilibrium become suspect. 

The results of two typical calculations are shown in figures 17 and 18. The differ- 
ence in the shape of the enthalpy distribution for the two cases is a result of the difference 
in nongray character of the radiation. In the low-temperature case (fig. 17, 

T s = 11 210° K), 95 percent of the radiant flux incident on the surface comes from 



(a) Enthalpy distribution. 

Figure 17.- Results for altitude = 54.85 km, 11.67 km/sec, 

RN = 300 cm, T s = 11 210° K, and p s /p 0 = 8.24 x 10"3. 
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optically thin portions of the spectrum. As a result the enthalpy distribution is similar 
to those for optically thin gray gases (see, for example, the curves for kp = 0 and 0.3 
in fig. 5(c)). However, at the higher temperature (fig. 18, T s = 14 021° K) only 45 per- 
cent of the radiation flux incident on the surface comes from optically thin portions of the 
spectrum. As a result the enthalpy distribution appears more like those for a gray gas 
with Bouguer number of order 1 (see curves for kp = 1.0 and 3.0 in fig. 5(c)). 

Spectral distributions of the radiant flux incident on the body surface are presented 
in figures 17(b) and 18(b) for the two cases. These distributions are, of course, based on 
the step-function model and do not accurately represent the detailed spectra for air. 
Radiation from the lines in the vacuum ultraviolet (X = 2000 A) is represented 



(b) Spectral distribution of radiant flux incident on the body surface. 
Figure 17.- Concluded. 
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schematically, as the absorption-coefficient model does not specify the number, locations, 
or widths of the individual lines. Despite these inadequacies the spectral distributions 
provide some physical understanding of the nature of radiant energy transport in a non- 
gray, nonisothermal slab of gas. For example, the flux in the optically thick portions of 
the spectrum originates in the cooler regions of the shock layer adjacent to the body sur- 
face and is characterized by the temperature in these cooler regions. The flux in the 
optically thin portions of the spectrum comes from the entire shock layer and is charac- 
terized by its average temperature. 

The variation with free-stream velocity, altitude, and body nose radius of the non- 
dimensional rate of radiant heat transfer to a nonreflecting stagnation surface is presented 
in figure 19. It can be seen that the radiant heat-transfer rate varies linearly with nose 
radius only for very small values of the heating rate. At higher heating levels the radi- 
ant heat-transfer rate increases more slowly with nose radius. If the nose radius were 
allowed to increase indefinitely, while altitude and velocity remained fixed, the radiant 
heat-transfer rate would reach a maximum and then decrease asymptotically toward zero. 
This type of behavior has been demonstrated by the gray-gas analyses of references 2 
and 12. 



(a) Enthalpy distribution. 

Figure 18.- Resuits for altitude = 67.10 km, W QO = 15.24 km/sec, 
R N = 30 cm, T s = 14 021° K, and p s /p Q = 2.135 x 10“3. 
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Because the radiative properties of high-temperature air are usually presented in 
the form of global (wavelength-integrated) fluxes emerging from one side of an isothermal 
slab of finite thickness, it is desirable to define a radiation cooling factor. This factor 
F c is the ratio of q R , the radiant heat-transfer rate for a self-absorbing nonisother- 
mal shock layer, to the heat-transfer rate for a self -absorbing isothermal, or adiabatic, 
shock layer. The influence of self-absorption nearly cancels in this ratio so that F c 
depends primarily on the fraction of energy lost by radiation from an isothermal shock 
layer and the enthalpy dependence of the absorption coefficient. The fraction of energy 
lost by radiation from an isothermal shock layer is proportional to a , the nondi- 

mensional global flux emerging from one side of an isothermal slab with thickness equal 
to the shock standoff distance. 



lb) Spectral distribution of radiant flux incident on the body surface. 
Figure 18.- Concluded. 
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(a) Altitude = 30.48 km. 

Figure 19.- Variation of the stagnation-point radiant heat-transfer rate with altitude, 
velocity, and nose radius. 
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(b) Altitude = 36.58 km. 


Figure 19.- Continued. 








<0 Altitude = 60.96 km. 
Figure 19.- Continued. 
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(h) Altitude = 73.15 km. 
Figure 19.- Concluded. 
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Since the absorption-coefficient model does not vary drastically over the range of 
conditions of interest, ti. enthalpy dependence of the absorption coefficient can be char- 
acterized by (though certainly not specified by) the enthalpy dependence of the Planck 
mean absorption coefficient at normal-shock equilibrium conditions (dtcpj dh) This 

quantity was found to be relatively constant over the range of altitudes for a given velocity 
except at 7.92 and 9.14 km/sec. plot of (dKpj dhj against velocity for various 

altitudes is presented in fig. 20.j 

Because velocity is a more convenient quantity than ^d/Cp^dhj^ the cooling fac- 
tors are shown as plots of F c against <lp w a for various velocities (fig. 21). It can 

be seen that the results correlate well on this basis except for the two lowest velocities, 
7.92 and 9.14 km/sec. The results from two other nongray studies (refs. 3 and 13), which 
neglected line radiation, are indicated in these plots by solid symbols. In two cases it 



Altitude, km 

O 30.48 
□ 36.38 

O 42.67 
A 48.77 
k 34.86 
k 60.96 
O 67.06 
O 73.15 




10 11 12 13 14 lb 16 

W M , km/sec 


Figure 20.- The enthalpy variation of the absorption coefficient as a function of 
altitude and velocity. 
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was necessary to adjust the results to a different velocity. This was done by plotting the 
values of F c computed herein against velocity for constant qp w & . The value of F c 

for the desired velocity was then obtained from a curve drawn parallel to the data, through 
the point obtained from the reference. The excellent agreement between the points 
obtained from references 3 and 13 and the correlations of this paper seem to indicate that 
the neglect of line radiation has little effect on the parameter (d/cp^dh) 

The results presented herein also correlated reasonably well when the cooling fac- 
tor F c was plotted as a function of the radiation cooling parameter e for all the con- 
ditions studied. This correlation is indicated by the hatched band in figure 22. Also 
shown in this figure are the results of a number of gray -gas calculations which have been 
presented in the literature. The results from reference 14 are restricted to a transpar- 
ent gas and were computed numerically by using integral methods. The results from ref- 
erences 2 and 15 include small amounts of self-absorption. Advantage was taken of the 
small values of Bouguer number kp in both studies to obtain solutions up to a 
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(b) Woo = 9.14 km/sec. 

Figure 21.- Cooling-factor correlation. 
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Figure 21.- Continued. 
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Figure 22 - Variation of the cooling factor with the radiation cooling parameter. 

quadrature, which were then evaluated numerically. The results from the references 
agree well among themselves ^kp and /cp are the parameters which produce the 
spreading of the curves) and lend support to the nongray results of the present study. 
Since the gray and nongray results show the same trends for small values of e, they 
might be expected to show similar trends for large values of e. If this is the case the 
data of the references can be used as a guide for extrapolation of the nongray results. 

Effect of Radiation on Convective Heat Transfer 

Even though the analysis of this investigation is based on the assumption that the 
gas in the shock layer is inviscid and a nonconductor of heat, some conclusions can be 
drawn regarding the effect of radiation on convective heat transfer. It can be shown that 
to first order in the boundary -layer parameter (Np e ^~*/2 the convective heating rate 

is proportional to the driving enthalpy difference across the boundary layer. (See, for 
example, appendix A of ref. 2 , which holds for the case of an emitting and absorbing gas 
because the enthalpy is an analytic function of the Dorodnitsyn variable 77 even at the 
wall, and ref. 16 , which holds for an emitting and nonabsorbing gas with nonanalytic 
enthalpy in the neighborhood of the wall.) If the wall is cold (as has been assumed 
throughout this investigation) the enthalpy of the wall can be neglected and the convective 
heating rate becomes proportional to the enthalpy at the outer edge of the boundary layer. 
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The location of the outer edge depends on the boundary -layer parameter and is given 
approximately by 



Values of the enthalpy h e at the edge of the heat-conducting boundary layer were 
obtained for the range of altitudes, velocities, and nose radii for values of the Peclet num- 
ber Npe of 10^, 1()3, and 10^. The results are presented in figure 23 as a function of 
the nondimensional adiabatic radiant heating rate for various velocities. A good correla- 
tion is obtained for h e for the same reasons that a good correlation was obtained for the 
cooling factor. 

The quantity h e is an approximation to the ratio of convective heating rate for a 
radiating shock layer to that for a nonradiating shock layer. When radiant energy 
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Figure 23.- Correlation of the boundary-layer edge enthalpy. 
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transport is important, the convective heating is reduced from the radiationless value 
(he = 1). If the analysis could be continued to indefinitely large values of the nondimen- 
sional adiabatic radiant-heating rate, the value of h e would become asymptotic to zero. 
This variation has been indicated by the gray-gas analyses of references 2, 12, and 17. 

Of course, figure 23 gives only an order-of-magnitude estimate of the radiation- 
convection coupling. Not included are the effects of variable transport properties, 
enthalpy gradients at the edge of the boundary layer (a second-order effect), and differ- 
ences in the characteristic Reynolds and Prandtl numbers between radiating and nonradi- 
ating cases. Also no account has been taken of radiation in the boundary layer. In the 
cool region of the boundary layer near the wall the gas will absorb more radiant energy 
than it will emit. This will tend to increase the slope of the enthalpy distribution adjacent 
to the wall, which, in turn, will increase the convective heat transfer somewhat. 

While absorption in the boundary layer will increase the stagnation- point convective 
heating rate it need not increase the stagnation-point total heating rate. It is apparent 
that a radiating gas in the boundary layer will (1) reduce the rate of radiant heat transfer 


h 

e 


i.o r 


o o m 


o io| 

□ 1CT 

o io 4 


.6 ' 
10 


-4 



(c) Woo = 10-67 km/sec. 


1.0 | 


.8 h 


.6 


« # # 


o 




°5V % 

v B 


o 


.4 


10 ' 


-4 


10 ' 


-3 


10 ' 


-2 


10 ' 


-1 


^R,w,a 


(d) Woo = 12.19 km/sec. 
Figure 23.- Continued. 
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to the body surface because of absorption, (2) increase the rate of convective heat transfer 
to the body surface because of the higher enthalpy level and resulting steeper enthalpy 
gradient at the surface, and (3) increase slightly the amount of energy convected out of 
the stagnation region by the tangential flow. If the effect of radiation in the boundary 
layer on the external inviscid flow is ignored (for conditions of interest here, it is cer- 
tainly a higher order effect) , the energy influx from the external inviscid flow to the 
boundary layer will not depend on the radiating character of the boundary-layer gas. 
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(f) Woo = 15.24 km/sec. 
Figure 23.- Concluded. 



Whether or not ihe increase in convective heating will be offset by the decrease in 
radiant heating depends . ' a large extent on the surface properties of the body. If the 
body surface were perfectly reflecting, the total heating rate to the surface would be 
increased by absorption in the boundary layer because only the convective component con- 
tributes to the total heating rate. If, on the other hand, the wall were perfectly ab sorbing , 
then consideration of an energy balance across the boundaryJLayer would show that the 
total heat transferred to the surface by radiation and convection will be less for an 
absorbing boundary layer than for a nonabsor bing b oundary layer because of the slightly 
larger amount of energy convected out ofjhe stagnation regi on by the tangential flow in 
the former case . 

For gases, most of the absorption by the boundary layer would occur in the vacuum - 
ultraviolet portion of the spectrum where, for most materials, the surface reflectivity is 
small. Thus, one might expect that the combined radiation and convection heat transfer 
to the surface as predicted by the inviscid radiating analysis and the nonradiating- 
boundary- layer analysis would be in reasonable agreement with, or perhaps slightly 
higher than, what would be predicted from solutions of a radiating boundary layer properly 
matched to the external inviscid solution. 

Radiation Blockage Effect of Ablated Vapor Layer 

When the radiant heat-transfer rates are large the ablated vapor layer will ab sorb 
more radiant energy than it emits because the characteristic shock-layer tempe rature 
will be considerably higher than the characteristic ablated-va por-layer _temp eratur e. If 
this is the case, the hot-air shock layer can be solved independently of the ablated vapor 
layer (such as has been done in this study) and the results can be used to establish bound- 
ary conditions at the outer edge of the ablated vapor layer. Solution of the ablated- vapor - 
layer problem is complicated by the importance of viscosity, thermal conductivity, species 
diffusion, and chemical reactions, and no attempt has been made to obtain a solution here. 
However, it has been pointed out that the ablated vapor layer will act to inhibit the trans- 
fer of radiant energy from the hot-air shock layer to the body surface. Since this is a 
benefi c ial e ffect, it is_ advisable (strictl y from th e point of view of redu cing the ra diant 
heat-transfer rate) to provide for an ablated vapor with large and ext ensive (in. wave length ) 
abso rption cross se ctions s uc h as those of the alkali metals (lithium, sodium, and 
potassium). 


Precursor Effect 

Some of the radiation energy emerging from the front of the shock layer is absorbed 
by the cold free-stream air in the path of the vehicle. Much of this absorbed energy 
ionizes nitrogen and oxygen molecules and dissociates oxygen molecules (some of the 
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resulting oxygen atoms may in turn be ionized by the strong vacuum -ultraviolet radiation). 
Because of the low density of the free-stream air it seems likely that the number of ions 
and electrons carried into the shock layer would be considerably larger than would be 
predicted for thermodynamic and chemical equilibrium. For the same reason it would be 
likely that the heavy particles (molecules, atoms, and ions) entering the shock layer would 
have very little thermal energy (characterized by a temperature not much in excess of the 
ambient temperature), while the electrons would carry more thermal energy (character- 
ized by a somewhat elevated electron temperature). The energy absorbed in the free 
stream would only slightly affect the density and velocity. Thus the energy absorbed in 
the path of the vehicle is returned to the shock layer largely in the form of chemical 
energy and thermal energy of electrons. Under the assumption of local thermodynamic 
and chemical equilibrium within the shock layer, the return of this energy will result in a 
value of enthalpy behind the shock larger than that predicted for no free-stream absorp- 
tion. It should be noted that the presence of electrons (even trace amounts) in the free 
stream will speed up the ionization processes behind the shock and increase the range of 
validity of the assumption of equilibrium in the shock layer. 

The momentum and energy balances across the shock wave are 

P s = \ P„W« 2 (40a) 

h s*'| w « 2 (l+yq R , s ) <40b) 

where p s is the pressure immediately behind the shock, h s * is the enthalpy immedi- 
ately behind the shock when absorption by the free stream is taken into account, g 
is the nondimensional radiant energy flux emerging from the shock layer in the upstream 
direction, and y is the fraction of this flux which is absorbed by the cold free-stream 
air and returned to the shock layer. The energy balance was obtained under the reason- 
able assumption that the free-stream air absorbs but does not emit radiant energy. The 
dominant effect of absorption by the free stream is an increase in shock-layer enthalpy, 
with no change in shock-layer pressure. Therefore the effect on the radiant and convec- 
tive stagnation-point heating rates can be estimated from results that do not include free- 
stream absorption by using an effective free-stream velocity Woo* and an effective free- 
stream density : 




* 


P 


oo 


s woo \f 1 + rq R,s 

(41a) 

*.pj w :) 2 

(41b) 

\WcoV 
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When q^ g is small, as it must be in the small-perturbation analysis, it can be con- 
servatively approximated by w a , the adiabatic heating rate. Also because of the 

smallness of q^ 0 , expressions (41) become 
r ,w,a 

W„* * W„ (l + \ yq R w (42a) 

< 42b > 

If Woo - Woo and p^ - p^ are small compared with W^ and p w , respectively, 

the rate of radiant heat transfer to the stagnation point with free-stream absorption 
becomes 



By combining equations (42) and (43) and ignoring terms of order ^yq-j^ w a) ^ or higher, 

a factor for the increase of radiant heat-transfer rate with free-stream absorption can be 
defined. This factor is 
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(44) 


This factor with y = 1 has been computed for the range of altitudes, velocities, and nose 
radii considered in this study, and the results are presented as functions of w a f° r 

various velocities in figure 24. These values of 3?-^ represent upper bounds because y 
will be less than 1. The determination of y is complicated by nonequilibrium thermo- 
dynamics and chemistry in the absorbing free stream and has not been attempted herein. 
Values of 3?-^ for y different from 1 can be obtained from the plots of figure 24 by 
means of the formula 


A factor for the increase in convective heat-transfer rate with free-stream absorp- 
tion can also be defined: 



The term on the far right-hand side of equation (46) utilizes the usual correlation of 
stagnation -point convective heat-transfer rates in terms of powers of the free-stream 
density and the free-stream velocity. Hoshizaki (ref. 18) gives values for the exponents 
of s = 0.5 and t = 3.2. These values are used herein. For Woo* - W M « Woo and 

™R,w,a <<: 
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(b) = 9.14 km/sec. 

Figure 24.- Correlation of the factor for increase in radiant heat transfer with free-stream absorption. 
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Thus 


„ „ i (d In h e \ 

~ 1 + ^R^a • + 2 \d i n WoJ p 


8 In h e 

9 In p~ 


On physical grounds it can be argued that h e will not increase with W TO and that 

- Wco 2 h e will not decrease with W^. The results of this investigation substantiate this 
2 

argument except for some cases where q„ is large and the expected precision for 

^ rijW jet 

h e is poor. Thus the first derivative in formula (47) is bounded as follows: 


^8 In h e \ < 

V 8 ln Woo/p " 
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(f) W M = 15.24 km/sec. 
Figure 24.- Concluded. 




Figure 25.- Bounding curves of the factor for increase in convective heating 
with free-stream absorption. 


Examination of the results indicates that the other derivative is bounded as follows: 


< 9 ln h e 

\9 ln Par 


Consequently, the factor 3>c is also bounded: 


1 + 0-1 ^R,w,a = # C ^ 1 + 2 •^R.w.a (48) 

The two limit curves for <3>p are shown as functions of yq^ in figure 25. Also 

ii, w, a 

shown are values of obtained in this study for y - 0.25, 0.50, and 1.00 for two dif- 

ferent combinations of free-stream velocity, altitude, and nose radius. 

Body Surface Reflectivity 

In the nongray studies the surface has been assumed to be a perfect absorber. This 
is a conservative assumption which was necessitated by the lack of reliable information 
concerning the radiative properties of high-temperature ablating materials. It is quite 
likely that a real surface will be partially reflecting in the visible and infrared portions of 
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the spectrum. The small-perturbation method presented in this paper can easily be 
modified to account for nongray surface properties by introducing a step function (in wave- 
length) model for the surface reflectivity. It is advantageous to specify the steps so that 
the edges coincide (in wavelength) with the edges of the step-function absorption- 
coefficient model. 

From the point of view of reducing radiant heat transfer to the body surface, a 
highly reflecting surface would be desirable. If the blockage effect of the ablated vapor 
is considered in addition to the radiative surface properties, it is apparent that an attrac- 
tive ablation material would be one which in the vapor form absorbs strongly in certain 
spectral regions (such as the vacuum ultraviolet), and in the solid, or liquid, form reflects 
strongly in the remainder of the spectrum. Whether or not a material exists which exhib- 
its a definite advantage over others is not now known. 

RESUME 

A singular perturbation solution to the blunt-body stagnation flow of an inviscid, 
radiating gas has been obtained by means of the P-L-K, or perturbation-of-coordinates, 
method. A number of results for a gray gas have been presented in order to provide 
some physical insight into the effects of various parameters on the shock-layer enthalpy 
profiles and the radiant heat-transfer rates. 

A nongray absorption- coefficient model was developed which includes, in an approx- 
imate way, the important vacuum -ultraviolet contributions of bound-free and line transi- 
tions. This model was used to obtain solutions pertinent to the case of reentry into the 
earth's atmosphere. While the results are restricted to small values of the radiation 
cooling parameter, which characterizes the relative importance of radiation and convec- 
tion as energy-transport mechanisms, they cover broad ranges of vehicle velocity, alti- 
tude, and nose radius, which are of practical interest. 

The characteristic enthalpy variation of the model absorption coefficient was found 
to be nearly independent of altitude and nose radius for fixed vehicle velocity except for 
velocities lower than 10.67 km/sec. Thus it was possible to correlate certain quantities 
by plotting these quantities as functions of the nondimensional adiabatic radiant heat- 
transfer rate for various altitudes and nose radii at fixed vehicle velocity. Among the 
quantities correlated was the cooling factor (the ratio of the stagnation-point radiant 
heat-transfer rate to the adiabatic radiant heat-transfer rate). The cooling-factor corre- 
lation is particularly useful because it eliminates the need to perform nonadiabatic cal- 
culations whenever radiant heat-transfer rates are desired. Also correlated was the 
factor by which the convective heat-transfer rate is reduced because of radiation losses 
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in the shock layer. Finally, upper bound estimates were made of the effects of absorption 
of precursor radiation by the free-stream air on the radiant and convective heat-transfer 
rates. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 24, 1968, 

129-01-02-03-23. 
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APPENDIX A 


MATHEMATICAL DEVELOPMENT 

In this appendix, the method of small perturbation will be used to obtain a solution 
to the integrodifferential system of equations governing the flow in the inviscid region of 
a radiating shock layer. Mathematical details that are not considered to be appropriate 
to the body of the paper are included herein. 

The system of equations to be treated is presented below: 


i(v)h'(v) + eq R '(r?) = 0 

(Al) 

Q 

2i(rj)i"(ri) - £f’(7j)J + a 2 h(? 7 ) = 0 

(A2) 

h (^ A ) = 1 

(A3) 

f (0) = 0 

(A4) 

f (^A) = 1 

(A5) 

f’(r? A ) = -j=i= 
\ 1 \ 2y(1 - y) 

(A6) 


Conventional Perturbation Procedure 

If the functions h(? 7 ,e) and f(??,e) are assumed to be analytic in the vicinity of 
e = 0, they may be written in the expanded form: 

OO 

h(? 7 ,e) = ^ e n h n (? 7 ) (A7) 

n=0 

oo 

e n f n (77) (A8) 

n=0 


It is anticipated that the first few terms of these expansions will provide an accurate 
estimate to the solution of the system of equations (Al) to (A6) when the parameter e is 
small compared with 1. 
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The integral term q^’ and the constant r also depend on the parameter e 
through their dependence on the function h (rj,e). These quantities are also assumed to be 
analytic functions of e near e = 0, so that 

oO 

= X eT ^R,n {7l) (A9) 

11=0 

OO 

^A= l e ^A,n (A10) 

n=l 


Substituting the expansions (A7) to (A10) into system (Al) to (A6) yields 
+ e fjho’ + ifli + q Rj o'(T7) + e 2 f2 h o’ + f l h l' + f 0h2’ + Vr,!^) + • • • = 0 (All) 



2f 0 f 0 " " (*0 + a2h 0 + 6 2*0*1 “ 2f 0 *1 + 2f 0”fl + a2 ^i 

+ e 2 2f 0 f 2 ” - 2f 0 'f 2 " + 2f 0 ”f 2 + 2 ^ 1 " - fix') 2 + a 2 h 2 + . . . = 0 (A12) 



63 


APPENDIX A 


Since the small parameter e is arbitrary, system (All) to (A16) can be satisfied 
only if each coefficient of each expansion in e is identically zero. This leads to a 
recursive set of purely differential systems. 

The zero-order system is 


5T 

II 

O 

(A17) 

2f 0 fo” - (f 0 ’) 2 + a2h 0 = 0 

(A18) 

h o(^A,o) = 1 

(A19) 

f 0 (0) = 0 

(A20) 

i— 1 

II 

<f 

o 

(A21) 

0 ( r7A >°) px(l - x) 

(A22) 

The solutions to this system are easily found. They are 


h 0 = 1 


f 0 = (1 - a )t] 2 + ar? 

(A23) 

^,0 = 1 

(A24) 

2\f2xd - X) 
a = 1 

(A25) 

1 + V 2 xd - x) 



The systems of first and second order may be written in the general form 


f 0 h n' + f l h n-l +< lR,n-l _0 

(A26) 

fQ^n" “ ^0 *n + ^0 *n = ^n^ 7 ?) 

(A27) 
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h n(!) = 4> n 

(A28) 

fn(0) = 0 

(A29) 

f n(7A,o) = ° n 

(A30) 

: n’K,o) = 

(A3 1) 


Equation (A26) can be integrated directly to obtain 


„ r lq R,n-l ,(x) +f 1 (x)h n . 1 *(x) 

h n(v) = 4>n + \ 2 T-7T 

J v f o( x ) 


the homogeneous solutions of equation (A27) are 77 + — and 77 2. As shown by 

2(1 - a) 

Ince (ref. 19) , the complete solution is 


^ - i [ 2(1 . a), ♦ - \ 

[(1 - a)x + a] 77 x 2 [(l - a)x + a] 


<I>n(x)dx 


2 Jo [d-a)x + a ] 2 


(A3 3) 


Substituting this expression into condition (A31) provides the following relation for 
the determination of 77 ^ : 


J o 


1 $ n (x)dx 


0 [(1 - a)x + a] 

The quantities <p n , and are 

®l(v) = a2hi(77) 


(A3 4) 


*2fa) = ~ f l f l" + |( f l') 2 " \ a2h 2 0?) 
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01 = 0 

^2 = 


6»! = -(2 - a)?7 A j 


e 2 .-L + ?^ 

I v 


1 $j(x)dx 


< 

a)x + aj ' 


^A,l " (2 " a)rj A,2 " (1 - a ^A,l 2 


01 = -2(1 - a)r7 A j 


1 ^(xjdx 


0 2 = -<2(2 - a)r/ A 1 + (2 - a)f — 

\ ’ J ° [d 


a)x + aj 


2)^,1 " 2 < J - a ^ A ,2 


(A3 7) 
(A3 8) 

(A3 9) 
(A40) 
(A41) 
(A42) 


The divergence of the radiant flux vector q R ’(^?) ma y be written in expanded form 

by substitution of (A7) to (A10) and the expansions of the quantities /c x (h), B x (h), 

T^(ri,e), and E n |V x (?7,e) into expression (7) of the body of the paper. For completeness 
the expanded forms of k x , B x , t x , and E n ^V x ^ are given here: 


+ • . . 


K x {h) = K X { 1) + €K x (l)hi(?7) + e 2 j^K x (l)h2(r?) + \ « x (l)hl 2 (»7) 

B x (h) = B x (l) + eB x (l)hi(?7) + e 2 j6 x (l)h2(rj) + 1 B x (l)h 1 2 (?]) + . . . 


rv 

V (v,e) s J o K x (7])d7] 

= k x (1)v + e/t x (l)j^ h l(v)dv + ■ 


(A43) 

(A44) 


(A45) 
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E n kpk x (r?) - r x (g) = E n /kp 


k x (1)(v - I) + efc x( 1 )^ 71 h 1 (x)dx 


= E n [k P K x (l)(7j - £)] 
- ekp/c x (l) 


n 7] 

, r 

hj(x)dx 

E n-ljkpKx(l)(7? - £,) 


(A46) 


The following property of the exponential integral function was used to obtain (A46): 

E n -l( x ) = E n (x) 

With these expansions in hand, expressions for the terms q\v) and Qp ^ ’ ( 77 ) 
can be obtained. The results are 


= "j* 0 K A B x[ E 2[ k x( 1 " *)] + ( X " r 0,x) E 2( k A^ dX (A47) 

= ~ hl ^ 0 ^X B ^|E2[kx(l - V)] + ( x - r O,x) E 2( k X 7 7)j’ “ 

+ £ h l(^ E l( k xb - *|) d « + ^X^E^d - r?)]^ h^d* 

+ ( x “ r 0,x) E lMlJ h l^ d *} + r l,X E 2 dx ( A48 ) 

In these expressions, the notation has been simplified somewhat by omitting the 
argument 1 in the terms k^, k B^ ? and and by introducing the quantities 


k A = k P*A 


(A49) 


r 


0,A 




(A50) 
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r l,A = 2r w 


«***£ h 1 (|)E 2 (k x ?)d? + K X B x E 2 (k x ) h!(«)d« 


+ j 7a,1 k X B A E 2 



P-L-K Solution 

It has been pointed out in the body of the paper that the first-order solution for the 
enthalpy distribution has a logarithmic singularity at the point 77 = 0 and the second- 
order solution behaves like the logarithm squared. As a consequence, the assumed 
expansion diverges as the origin is approached and the small-perturbation solution is not 
uniformly valid. In order to obtain a solution which is uniformly valid throughout the 
domain of the problem, the Poincare-Lighthill-Kuo method (see ref. 7) is used. In this 
method the independent variable as well as the dependent variables is expanded in a 
McLauren series of e. For this problem, 

7 j = x + e 77 j*(x) + e 2 7 j 2 *(x) + . . . (A52) 

h( 7 ?,e) = h 0 *(x) + eh!*(x) + e 2 h 2 *(x) (A53) 


The asterisk is used here to distinguish between the coefficients in the P-L-K expansion 
and the coefficients in the conventional expansion (eq. (A7)). The quantities f ( 77 ) and 
q^'(? 7 ) may also be expanded in terms of x as follows: 

00 

f (v) = ^ e n f n (77) 

n =0 


S X enf n[ x + er ?i*( x ) + • • 

n =0 

= f()(x) + e[fi(x) + ^i*( x )fo'( x )] + e 2 jf 2 (x) + ^(xJfi'O?) 
+ J7 2 *(x)f 0 '(x) + i[?7 1 *(x)] f 0 "(x)J + . . . 


(A54) 


V (7?) =q R,o' (x) +eq R,l’ (x) 


^!*( x ) q R,0 ^ X ) 


(A55) 


When expansions (A52) to (A55) are substituted into equation (Al), a set of equa- 
tions for hQ*(x), h 1 *(x), h 2 *(x), . . . results. The quantities 7 j 1 *(x), ?? 2 *(x), . . . 
and their first derivatives also appear. These quantities are arbitrary and should be 
chosen in such a manner as to reduce the strength of the singularities in the terms 
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h 2*( x )> h 3*( x )> • • •> s that these singularities are never stronger than that of the low- 
est order singular term v ■ the present case, the first-order term). Pritulo (ref. 8) has 
shown that the coefficients in the expansion of h(ry,e) in the P-L-K method are related 
to the coefficients of the conventional expansion in the following manner: 

h 0 *(x) = h 0 (x) (A56) 

hl*( x ) = h^x) + 7? 1 *(x)h 0 , (x) = hj(x) (A57) 

h 2 *(x) = h 2 (x) + T 7 1 *(x)h 1 ’(x) + i? 2 *(x)h 0 '(x) + |[? 7 1 *(x)] 2 h 0 "(x) 

= h 2 (x) + Tjj^hj^x) (A58) 


Now, instead of choosing differential equations for the ?? n *(x) to satisfy the criterion 
previously stated, the values of the 7? n *(x) can be chosen directly. In this case, an 
obvious choice is simply 


7 7 1 *( x ) = - 


h 2 ( x ) 

hi ? (x) 


(A59) 


which gives h 2 *(x) = 0. ^This choice satisfies the condition ^*(1) = 77 ^ j.J 

Transforming the independent variable by means of formula (A52) removes the sin- 
gularity from the domain of the problem. That this is true can be seen by noting that the 
condition 77 = 0 does not imply x = 0 but rather (for this problem) implies that x is 
some small positive number. Hence, to first order, 11 ( 77 , e) = 1 + ehi*(?7) is nonsingular 
in the domain 0 § 77 = 1. 


Method of Matched Asymptotic Expansions 

V&n Dyke has pointed out (ref. 20) that the method of matched asymptotic expansions 
is applicable whenever the P-L-K method can be used. Thus, it would be interesting to 
formulate the solution when radiation is a small perturbation by using the method of 
matched asymptotic expansions. Use of this method implies that the domain of the prob- 
lem can be divided into at least two regions in which the governing equations take on dif- 
ferent asymptotic forms. There must also be some overlap between adjacent regions so 
that a smooth transition between solutions valid in these adjacent regions can be effected. 
In the problem of this chapter, the regions are the "outer” region in which the conven- 
tional perturbation solutions are valid and the "inner” region in the vicinity of the wall at 
77 = 0. 
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The equations which describe the conditions in the outer region are simply the sys- 
tem (Al) to (A6). In order to obtain the "boundary-layer" form of these equations, it is 
necessary to stretch the coordinate 77 in the vicinity of the wall. This stretching takes 
the nonlinear form 


no = [f(T?)] e 

F'(?) = f’(r?) 


(A60) 


where £ is the stretched boundary- layer coordinate and F(£) is the velocity function 
written in terms of £. It follows that 


= e (A61) 

*) f(t|) 

and the energy and momentum equations, respectively, take the forms 

F(?)H’(?) + q R '(?) = 0 (A62) 

2F(C)F"(?) - [f’( 0] + a2H(?) = 0 (A63) 

where 

H(0 = h(?7) (A64) 

and 

q R ’(?) = q R '(?7) ( A65 ) 

One boundary condition is available: 

F(0) = 0 (A66) 

The remaining two constants of integration can be obtained by matching the inner and outer 
solutions according to the matching principle set forth in reference 20. 

The boundary-layer system is seen to be quite complex. The energy and momentum 
equations remain coupled so that it is necessary to obtain a simultaneous solution to the 
two equations. Thus, as is often the case when the P-L-K method can be applied, its 
application is much simpler than the method of matched asymptotic expansions. 
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Step-Function Absorption- Coefficient Model 

If the spectral variation of the absorption coefficient is approximated by a step func- 
tion in wavelength, the wavelength integrals which appear in equations (A47) and (A48) can 
be replaced by finite summations with the following result: 


I 

'•r.o'*”' = - I - 1)] + (» - r 0,i) E 2( k i’')} 


(A67) 


q R,i' (7?) = " h l^) X " *?)] + (l - r 0,i) E 2(M)j> “ 


2 *i B i 


J kfUiBi ^ 1 hi($)Ei ki|r? - || d£ + KjBj^E i[ki(l - rj)] h^d* 


(l - r^E^ktf) h!(«)d«\ + r lji E 2 (k i 7 7 )j 


(A68) 


where jq is the value of the nondimensional mass absorption coefficient for the ith step, 
and 

Bi = B(\i_i,Xi) = f 1 B x dX (A69) 

v 1 J M-l 

Here X^_j and Xj are the wavelength values for the edges of the ith step. Also, 


kj = kp/q 


(A70) 


r 0,i r w,i 


r l,i = 2r w,i 


1 - 2E 3 (ki 

/CiBi^ h 1 («)E 2 (k l «)d| + /c i B i E 2 (k i )J o 1 hiU)d« + 77 A1 /c i B i E 2 (k i ) 


(A71) 

(A72) 
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MODEL NONGRAY ABSORPTION COEFFICIENT 

The various radiative processes of high-temperature air are arbitrarily divided 
into three groups herein for convenience of discussion. These groups are (1) atomic and 
ionic lines, (2) N + - and 0 + -electron recombination in the vacuum ultraviolet, and (3) others 
(which includes N + - and 0+-electron recombination in the near-ultraviolet and visible 
regions, free-free electron transitions, and molecular band systems). Only the processes 
of group 3 were considered in the computations of the radiative properties of air by Kivel 
and Bailey (1957, ref. 21), Armstrong et al. (1961, ref. 22), Riethof and Nardone (1962, 
ref. 23), and others. As a result of these calculations it was concluded that the shock 
layers of most entry vehicles would be optically thin. 

In 1963, Nardone et al. (ref. 24) and also Biberman and Norman (ref. 25) pointed out 
the importance of processes of group 2; that is, the N + - and 0 + -electron recombination 
in the vacuum ultraviolet. It was shown that this contribution is even more important 
than that of group 3 and, furthermore, the group 2 contribution invalidates the conclusion 
that most shock layers of interest are optically thin. In 1965, Hahne (ref. 9) and, inde- 
pendently, Sherman and Kulander (ref. 26) performed careful calculations of these group 2 
contributions. Their separate calculations showed good agreement. In 1964, Biberman 
et al. (ref. 27) showed the importance of the contribution of atomic and ionic lines 
(group 1) to the radiative properties of air. Allen (1965, ref. 10) included the line con- 
tributions for nitrogen and oxygen atoms (but neglected ionic lines, which can be impor- 
tant for combinations of high temperature and low density3) in his calculations of the 
radiative properties of air. 

For purposes of calculation in this paper, the absorption coefficient was obtained 
from what is considered by this investigator to be the best current information on the 
radiative properties of air. A description of the methods by which, and the sources from 
which, the model nongray absorption coefficient was obtained is given in this appendix. 

^Vorobyov and Norman (ref. 28) show that the lines of the N+ ion in a nitrogen 
plasma begin to become important at temperatures as low as 8000 ° K for a density ratio 
p s p/Po = 10-3. For a density ratio PsP/P 0 = 10-2 the N + ion lines begin to become 
important for temperatures greater than about 12 000° K. It should also be noted that the 
f-numbers used by Vorobyov and Norman in their calculations for lines in the ultraviolet 
are, in many cases, an order of magnitude larger than those used by Allen (ref. 10). The 
scarcity of experimental evidence, combined with the large uncertainties with regard to 
line width calculations, makes it difficult to assess which set of f-numbers is the more 
reliable. The results of Allen have been used for the calculations of this paper because 
they are more extensive than those of Vorobyov and Norman. 
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Line Radiation 

The contribution of the atomic lines to the absorption coefficient of air was obtained 
from the data of Allen (ref. 10). The neglect of the ion lines by Allen is not believed to be 
a serious source of error for the calculations of this paper. Allen presented the data on 
line contribution in the form of spectrally integrated radiation fluxes emitted from one 
side of isothermal, infinite slabs of finite thickness (specifically 0.1, 1, 5, 10, 15, and 
20 cm) for all lines with wavelengths greater than 2000 A. Presentation in this manner 
obscures the spectral characteristics and makes it difficult to obtain the monochromatic 
absorption coefficient. 

A study of Allen's data indicates that all lines with wavelength greater than 2000 A 
are optically thin for values of temperature, density, and slab thickness of interest. As 
a result of this fortuitous circumstance, a Planck-type average absorption coefficient for 
these lines can be used in the calculations of this paper without any sacrifice in accuracy. 
(A Planck average is the proper wavelength average if, and only if, the gas volume is 
optically thin at all wavelengths.) The Planck-type average for lines of wavelength 
greater than 2000 A can be obtained from Allen's data by means of the formula 


Pk( 2000, 100 000) = 1 1( 2 000; *0° 000 ;> L ) per centimeter (Bl) 

oT s 4 L B(2000, 100 000) 


where I^\i,A 2 ,l) is the radiant flux for lines with wavelength between Aj and A 2 
emitted from one side of an infinite slab of thickness L (in centimeters). The quantity 
B(Ai,A 2 ^ is defined by the formula 

•*^2 


B 


(m,x 2 ) = f 


M 


B X dA 


(B2) 


Not all the lines for wavelengths less than 2000 A are optically thin for conditions of 
concern. Unfortunately, the manner in which Allen presents the information makes it 
impossible to determine the spectral characteristics of the line radiation; therefore, the 
stratagem described below is used to obtain a model for this contribution. First of all, 
it should be noted that the lines under discussion are restricted to the wavelength interval 
of 911 to 1800 A. (See ref. 11 for location of lines for the nitrogen and oxygen atoms.) It 
is assumed that the absorption coefficient is composed of a number of narrow steps (in 
wavelength) of uniform height superimposed upon a gray background. The widths and 
locations of the steps are arbitrary as long as the widths and the spacings between steps 
are smaller than the wavelength interval required to achieve a significant change in the 
value of the Planck function B^ and if the steps are distributed statistically throughout 
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the wavelength interval of 911 to 1800 A. The value of in the gray background is 
denoted by and that in the steps by o^. (See fig. 26.) 

When this model is used to compute the radiant flux emerging from one side of an 
infinite, isothermal slab of thickness L, the following formula is obtained: 


1(911, 1800, L) = aT 4 < |l - 2E 3 (pa 1 L) SB(AX background ) + 1 - 2E 3 (pa 2 E) 


SB (AX 


steps) 


(B3) 


With the previously specified restrictions on the widths, spacings, and distribution of the 
steps, this expression can be approximated by 


C r- 


1(911, 1800, L) = crT 4 B(911, 1800)^1 - 2E 3 (po' 1 l) 


+ (1 - 0)1 


- 2E 3 (pa 2 L)]| 


(B4) 


J J 



Here 1-/3 is a relative density for the steps. 
The three arbitrary parameters a\, cq, and /3 
introduced by this model can be evaluated by 
solving equation (B4) for three different values of 
slab thickness L, using Allen's data for the cor- 
responding values of 1(911, 1800, L). However, 
when this procedure was followed the resulting 
values for the parameters c^, a 2 , and /3 dis- 
played some anomalies which were believed to be 
associated with (1) the inability to read accurate 
values from the graphs in reference 10 and (2) 
the neglect, in that reference, of the effects of 
line overlapping on self-absorption. 


Figure 26.- Model absorption coefficient In order to avoid these difficulties, a pro- 

for ultraviolet lines. gram was developed in which the monochromatic 

absorption coefficient is computed for the nitrogen 
and oxygen lines in the ultraviolet. The line parameters (f-numbers, energy levels, and 
so forth) are the same as those used by Allen (i.e., listed in Griem, ref. 11). The expres- 
sions used to compute the line widths in the program were: 


For J(J + 1) ^ 


13. 6Z 2 


®ion “ ^upper 


w = 1.14 x 10 16 — iMZ? /_ 27.2Z2 


\[T Ei 


E, 


ion ^upper 


^ion " ®upper 


+ i 


(B5a) 
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For J(J + 1) < 


13. 6Z 2 


E 


ion 


E, 


upper 


w 


= 1.14 x 10 6 


N e 13. 6Z 2 

68Z 2 

\|T Ei on - E U pp er 

Eion " E U pper 


+ 1 - 3J(J + 1) 


(B5b) 


Here, 

J angular-momentum quantum number of the upper level of the atom for 

absorption 


Z charge of the atom plus 1 in units of the charge of an electron 

E ion ionization energy of the atom, cm“l 

E U pper energy of the upper level, or final state, of the atom for absorption, cm _ l 
w half of the line half-width, cm~l 

N e electron density, particles/cm^ 


Although sizable differences in the half-widths of some lines as calculated by this for- 
mula (B5) and the formula used by Allen were noted, these differences were essentially 
canceled when the individual line contributions were summed over all lines at wavelengths 
less than 2000 A. 

This program was used to compute the radiant flux emerging from isothermal slabs 
of various thickness L where only emission and absorption by the ultraviolet lines were 
considered. Results obtained from this calculation are presented in figure 27. The 
parameters oq, a< 2 ? anc * /3 displayed a great deal more regularity with variations of 
temperature and density when evaluated with the fluxes computed here than when they 
were evaluated with Allen* s data. 

It is expected that the step-function model for the absorption coefficient for the 
ultraviolet lines will give a reasonable estimate of the effects of absorption provided that 
the values of L chosen do not differ too greatly from the shock standoff distance and 
that the temperature distribution in the shock layer does not differ too much from 
isothermal. 
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N + - and 0 + -Electron Recombination in the Vacuum Ultraviolet 

The contribution to the optical properties of air of the N + - and 0 + -electron recom- 
bination in the vacuum ultraviolet was obtained from the data of Hahne (ref. 9). He indi- 
cated the existence of 12 thresholds, or photoionization edges. However, only four of 
them appear to be important. Consequently, the contribution to the absorption coefficient 
was approximated by four steps in wavelength with edges corresponding to the threshold 
wavelengths of 852, 911, 1020, and 1130 A. The step corresponding to the threshold at 
852 A was arbitrarily cut off at 400 A because the emission from air appears to be negli- 
gible for wavelengths shorter than this. The heights of the steps were set equal to a 
Planck-type average of the absorption coefficient over the interval encompassed by the 
step. The step-function model ignores the displacement of the thresholds due to inter- 
action effects (see ref. 27), which tends to shift the threshold to longer wavelengths with 
increasing ion concentration. 


Other Contributions 

The remaining contributions to the optical properties of air result from a combina- 
tion of free-free transitions, bound-free transitions, and molecular band systems. The 
principal of these, along with the wavelength interval in which each is important, are as 
follows: 


£(bound-free + free-free) 


>200 - 10 000 A 


Band systems 


O + e“ — O" + h^ 
N + + e~ — N + h^ 
0+ + e“ O + hi/ 

O free-free 
N free-free 
(N + + 0 + ) free-free 


N 2 (1-) 3000 - 6000 A 

N 2 Birge-Hopfield 1 and 2 900 — 2000 A 

N 2 (1+) 3000 - 20 000 A 

N 2 (2+) 2500 - 5000 A 


nJ Meinel 

5500 - 30 000 A 

NOj3 

2000 - 5800 A 

NOy 

1800 - 3200 A 

C >2 Schumann-Runge 

1800 - 5000 A 
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The maximum monochromatic absorption coefficient at wavelengths greater than 
1800 A is so small that the maximum monochromatic shock layer optical thickness is 
very much less than 1 for the temperatures, densities, and shock standoff distances of 
interest. Consequently, the contributions of the various radiators can be legitimately 
accounted for with a Planck-type average absorption coefficient for the wavelength inter- 
val of 1800 to 100 000 A. The contributions of the N 2 Birge-Hopfield systems were 
approximated with Planck-type averages over the wavelength interval indicated in the 
table. Although the shock layer for vehicles of interest is not optically thin in portions 
of this wavelength interval, the Planck-type averaging is not expected to introduce any 
significant errors into the analysis because absorption by these processes is small and 
relatively insensitive to wavelength in comparison with the group 2 contributions. 

Absorption-Coefficient Model 

Contributions from the three groups were combined to obtain a model for the 
absorption coefficient of air for a range of temperatures from 8000 to 16 000° K and 
densities from 10 - 3 to 10" * times the standard sea-level density. A compilation of the 
step-function model parameters for these conditions is presented in table I. An example 
of the model absorption coefficient is shown in figure 14. 
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